Orbital excitations in LaMnO^ 
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We study the recently observed orbital excitations, orbitons, and treat electron-electron correla- 
tions and lattice dynamics on equal footing. It is shown that the orbiton energy and dispersion 
are determined by both correlations and lattice-vibrations. The electron-phonon coupling causes 
satellite structures in the orbiton spectral function and the elementary excitations of the system are 
mixed modes with both orbital and phonon character. It is proposed that the satellite structures 
observed in recent Raman-scattering experiments on LaMnOa are actually orbiton derived satellites 
in the phonon spectral function, caused by the phonon-orbiton interaction. 
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Elementary properties of electrons in atoms and solids 
are determined by their charge and spin, but in many 
correlated electron materials also by their orbital degree 
of freedom ||l|,|[. In a free transition metal atom, with an 
open d shell, the 3d energy levels are five-fold orbitally 
degenerate, where each orbital state corresponds to a dif- 
ferent quadrupolar charge-distribution in real-space. In 
correlated Mott-insulators, among which there are many 
transition metal oxides (TMO's), monopolar charge ex- 
citations, that involve moving charge from one atom to 
another, are only possible at high energies because of the 
large Coulomb interaction between electrons. However, 
low energy multipolar charge excitations -corresponding 
to orbital excitations- are possible as they can be locally 
charge-neutral. 

In a solid the orbital degeneracy of a free ion is lifted. 
There are two physically distinct mechanisms. One pos- 
sibility is that the electron-electron interaction splits the 
states via the superexchange and in this way couples 
orbitals to the spin and relates orbital- to magnetic- 
order [pi. On the other hand, also the electron-phonon 
interaction, which gives rise to orbital order accompa- 
nied by a Jahn- Teller lattice distortion |^, can lift the 
degeneracy. The resulting orbital order is found in many 
TMO's, ranging from titanates (e.g. YTiOa Q), vana- 
dates (e.g. V2O3 il, LiVOa §, YVO3 [00 and mangan- 
ites (e.g. LaMnOa [|l, Ndi/sSri/jMnOs y), to cuprates 
(e.g. KCuFs 0). The quadrupolar charge ordering - 
orbital ordering- should give rise to elementary excita- 
tions with orbital signature, as the order causes a break- 
ing of symmetry in the orbital sector. The existence of 
such excitations, orbitons, was predicted in the 70's and 
theoretically studied ever since p|,p^-p2[ , but only very 
recently orbitons were observed for the first time in a 
Raman scattering experiment on LaMnOa [Ol. The ob- 
served orbitons were interpreted by some as being due 
to electron correlations ||l^ , but others suggest that they 
originate from the electron-lattice coupling WM. This 
motivates us to address the question of the origin of the 
orbitons, also because it is important to establish the 



precise nature of the orbitons as they, in turn, have a 
large effect on spin ^5 16 1 and monopolar charge excita- 
tions 0,|ll. 

We use a realistic model Hamiltonian for LaMnOa that 
incorporates both superexchange and electron-phonon 
coupling. We study the Hamiltonian first in the local- 
ized limit in order to gain more physical insight into the 
problem. This approach illustrates that for groundstate 
properties it is usually sufficient to treat the Jahn- Teller 
phonons as classical entities [|9|, but that for orbital 
dynamics it is essential to treat the lattice-vibrations 
quantum-mechanically. The full calculation shows that 
in LaMnOs the orbiton has exchange and lattice char- 
acter: its energy and dispersion are determined by both 
correlations and phonons. We propose that the peaks in 
the Raman-scattering data on LaMnOa |l^, are orbiton 
derived satellites in the phonon spectral function, which 
arise due to the mixing of the orbital and phonon modes. 

Hamiltonian. We consider the two-fold degenerate 
manganese eg states, with one electron per site. The 
electron can either be in the x'^ — y^ or 3z^ — r^ orbital, 
or in any linear combination of these two states. The 
interaction between neighboring orbitals is mediated by 
the superexchange and the electron-phonon (e-p) inter- 
action couples the electron to the two-fold degenerate 
Jahn- Teller phonons that have eg symmetry. We also 
take into account the dispersion of these phonons. Let 
us split up the Hamiltonian in an orbital, e-p and free 



+ rrQ2^, (1) 



where the sum is over neighboring sites (ij) along the 
r = a,b,c crystallographic axes. The orbital operators 
T^ can be expressed in terms of the Pauli matrices r: 
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rpa/b ^ ^^, _^ V3rf )/2, and T[ = rf. The e-p coupfing 
constant is denoted by g and the phonon operators of the 
so-called Q2 and Q3 Jahn- Teller modes with eg symmetry 
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where the local phonon energy is given by loq and the 
nearest neighbor coupling between the phonons by toi. 
The coupled Jahn- Teller modes along the three spatial 
axes are Q"/^ = {Q3, ± ViQ2^)/2 and Q', = Qg^ M- 
Note that the orbital excitations are locally charge neu- 
tral and do therefore not couple to breathing mode dis- 
tortions. 

The Hamiltonian Eqs.(y||) is general and one needs to 
make it more specific in order to describe LaMnOa. This 
is a prototype of an orbital ordered and Jahn- Teller dis- 
torted system and its crystallographic structure basically 
consists of corner connected MnOg octahedra, where the 
space between octahedra is filled with Lanthanum atoms. 
Due to the correlation and Jahn- Teller coupling the oc- 
tahedra are elongated, with the axis of elongation along 
the crystallographic a-direction on one sublattice, and 
along the 6-direction on the other, such that the recip- 
rocal lattice vector for the orbital order is Q = (tt, it, 0). 
An intersite phonon coupling arises because the elonga- 
tion of a MnOg octahedron induces an contraction of the 
neighboring octahedron, and vice versa, as the octahe- 
dra have one corner in common. The orbital order is 
formally incorporated in the Hamiltonian by performing 
a rotation of the T operators ||ll|] , with equal orbital ex- 
change constants Ja — Jb — J along the two axes in the 
plane. 

Transformations. In analogy with linear spin wave 
theory, the orbital modes can be found by performing 
a Holstein-Primakov transformation ||ll| . We introduce 
on each site i the bosonic orbital operators q[^ and qu: 
T- = \-qliqu and rf = \{q[i+qii)- We see from Eq. (|) 
that this transformation introduces a term that is lin- 
ear in the phonon mode Qs in H^-p- This is a conse- 
quence of the long range orbital order that we assumed 
to be present from the beginning: the lattice deforms ac- 
cording to the symmetry of the occupied orbital on each 
site. The linear term can be gauged away by introducing 
93 ~^ 93 +^7 7 where the shift rj is given by 77 = g/{L0Q—6uJi). 
After this shift we collect the quadratic and cubic terms 
in the bosonic operators and find in Fourier space 
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with Q^k = qlk + qv-k 

l\\k + l±kJc/J, l2k = 27±fe 
cosky)/2 and j±k — cosk^ 



and the dispersions 7n; 
- 7|H;, where 7||fc = {cosk^ 
PI. 
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FIG. 1. (a) Orbiton (Qi) and (b) phonon (Q2) spectral 
function in the localized limit. The first pole due to the orbital 
exchange is indicated at 3J. In the orbital spectrum the static 
crystal-field energy J is indicated by V. The spectral weight 
in (b) is multiplied by 10 for uj > ujo- 



Three important consequences of the orbiton-phonon 
coupling are present in Eqs. (Q^. First, the coupling to 
the lattice moves the orbiton to higher energy an amount 
4(777. This shift has a straightforward physically meaning: 
it is phonon contribution to the crystal-field splitting of 
the Cg states caused by the static Jahn- Teller lattice de- 
formation P2I. If, however, an orbital excitation is made, 
it strongly interacts with the Q3 phonon (Eq. H), so that 
the orbital excitation can be dynamically screened by the 
Jahn- Teller phonons and lowered in energy. The crystal- 
field splitting and screening are strongly competing as 
both are both governed by the energy scale set by the e-p 
coupling. Finally, the orbital and Q2 phonon modes mix, 
as is clear from the second term of H^-p- This implies 
that the true eigenmodes of the coupled orbital-phonon 
system have both orbital and phonon character. 

Localized limit. We illustrate the three physical effects 
of the orbiton-phonon coupling, discussed above, by con- 
sidering the Hamiltonian of Eqs. (|^H) first in the local- 
ized limit, neglecting all dispersion. The Hamiltonian 
then reduces to 

Hioc = {J + 2gQ3)qlqi + uJo{qlq3 + qlq2) + gQiQ2, (6) 

with J = 3 J + 4(777, which is the sum of the local or- 
bital exchange energy and static phonon contribution 
to the crystal-field splitting. The Hamiltonian with- 
out the last term is exactly solvable by a canonical 
transformation 123] so that we can obtain the expres- 
sions for the 6 by 6 matrix of bosonic Green's func- 
tions Du{v,t-t') = ~i{\Tql{t)q,{t')\), Dui^^t - t') = 
-i{\Tq^{t)q,{t')\), D2i[u,t-t') = -i{\Tql{t)ql{t')\), 
D22{i^,t-t') = -i{\Tq^{t)ql{t')\), withz/= 1,2,3. The 
last term in Eq. m) couples the orbiton and Q2 phonon 
mode and introduces a self energy pJ. We use the 
on-site Jahn- Teller vibrational energy uq = 80 meV 
as unit of energy, and J = uJo/2, in accordance with 
Refs. O.O]. In Fig. Ha the calculated orbiton spectral 
function, ——Im D(l,w), is plotted. For small e-p cou- 
pling g, most of the spectral weight is in the pole at 
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FIG. 2. Orbiton and phonon dispersion, neglecting dy- 
namical effects due to the e-p coupling; (a) without e-p cou- 
pling g and without bare phonon dispersion, (b) g/coo = 1/2, 
no bare phonon dispersion, and (c) g/coo — 1/2, finite bare 
phonon dispersion. The points of high symmetry in the Bril- 
louin zone correspond to those of Ref. [13] . 



Lu w 3J, and phonon satellites with decreasing inten- 
sity are present at higher frequencies, at energy intervals 
ujQ. The satellites are also known as Frank-Condon side- 
bands p^ , and their weight increases with increasing c-p 
coupling strength. For larger coupling constants the av- 
erage orbiton excitation energy increases, caused by the 
increase of the crystal- field splitting, but low and high en- 
ergy satellites are always present due to the interaction 
of the orbital excitation with lattice vibrations. 

The mixing of orbital and phonon mode gives rise to 
one extra phonon satellite in the orbiton spectral func- 
tion, at frequencies below ojq. In Fig. nb we see that, vice 
versa, due to the mixing a low intensity orbital satellite 
at w 3J is present in the Q2 phonon spectral function. 
The Q2 vibrational mode softens with increasing g, in 
contrast to the the Q3 mode, which is not affected by the 
e-p interaction. 

Full Hamiltonian. The quadratic parts of the Hamil- 
tonian in Eqs.(y-^ can be diagonalized by a generalized 
Boguliobov transformation. For the moment we do not 
consider the cubic terms, that give rise to the dynamical 
screening of the orbital excitation. In the transformed 
operators the Hamiltonian is Hguad = J^uk^t^kolka^k, 

with ql^ = E^w^M^Ife + ^^A'^M-fc- The energies of the 
eigenstates are given by 



'1,2k 



^3k 



Ck + xk± [(Cfe-xfe)' + 45^7^0] 

u;o{ujo - 67||A;CJi), 



(7) 



with Cfe = J{J - 7ife)/2 and Xk = ^0(1:^0 + 272fcWi)/2. 
The analytical expressions for uj;„ and w^„ are rather in- 
volved. In Fig. p^ the dispersion of the eigenmodes eik 
and e2fc are shown for different sets of parameters, where 
we used Jc/J — 0.5 |13|. The modes with predominantly 
phonon character are centered around E ^ loq, and with 
orbital character around the crystal-field energy J. For 
the parameters indicated in of Fig. 0b the wavefunctions 




V j b 
V 
- 1 lg=0.35 


- 1 1 ^' 


I *10 



12 3 

co/a)„ 



FIG. 3. (a) Orbiton spectral function at the F and 

X-point, g — (jJo/2. (b) Spectrum of the Raman-active Ag 
and Big phonon modes for loi/loo — 0.05 and g/uJo ~ 0.35. 
The experimental peak positions are indicated by V. For 
u! > ujQ the spectral weight is multiplied by 10, g = 0.35a;o, 
uji/ujo = 0.05. 



of the low energy excitations have on average 95% phonon 
character and the high energy modes 95% orbital char- 
acter. Note that the orbiton dispersion is almost entirely 
due to the exchange coupling J and that via the e-p cou- 
pling the orbiton dispersion reflects itself in the effective 
phonon dispersion. 

We now consider the interaction between the eigen- 
modes via the cubic term in Eq.(0). The cubic terms can 
be taken into account in a diagrammatic expansion. The 
first non-zero diagram corresponds to an orbiton that 
excites a phonon, propagates and absorbs the phonon 
again, which is a second order process. We calculate 
the self-energy due to this process self-consistently, i.e. 
instead of the bare orbiton propagator, we use the or- 
biton propagator dressed with phonon excitations. This 
is equivalent to the self-consistent Born approximation, 
as is used, for instance, to calculate properties of a sin- 
gle hole in the t-J like models ||l^,^. This approxi- 
mation works exceptionally well in t-J like models, for 
the problem of a single hole that is strongly coupled to 
magnons, or phonons |26[| , giving us confidence in its ac- 
curacy for the problem of a single orbital excitation cou- 
pled to phonons, which we consider here. 

After the Boguliobov transformation, the single cubic 
term in Eq.(0) maps into eight different non-zero cubic 
combinations of the new ai,k operators. We can make use 
of the observation that, in the parameter regime that we 
consider, the eigenmodes ei.2k have almost entirely or- 
biton c.q. phonon character. The calculated self-energy 
is therefore dominated by only one of the eight cubic com- 
binations of a^fe operators, which is at least an order of 
magnitude larger than other terms. Physically, this term 
is due to the orbiton-phonon scattering process described 
in the previous paragraph. Details of the derivation will 
be published elsewhere [P4| . 

The self-consistent calculation is performed numeri- 
cally, taking about 10^ points in the Brillouin-zone and 
an energy-grid with meshsize luq/IOO. The resulting or- 



biton spectral function at two high symmetry points in 
the Brillouin zone is shown in Fig. Ha, where we used 
the same parameters as in Fig. gc. Comparing these 
two figures, we see that due to the dynamical e-p cou- 
pling the poles with mainly orbital character are shifted 
to lower frequency and that at higher frequency phonon 
induced satellites develop. This is not unexpected, as the 
same happens for the system in the localized limit (see 
Fig. ^a). A closer look to the orbital spectrum at the 
F-point, however, reveals also that the effective orbital 
dispersion is wo/2, whereas the free orbital dispersion is 
luq- The reduced dispersion can be understood as a con- 
sequence of polaronic band-narrowing: the effective mass 
of the orbital excitation increases because of its dressing 
with phonons. 

Finally we can compare the calculated spectral func- 
tion of the Raman-active Ag and Big phonon modes, 
shown in Fig. |3|b, with experiment |13| . The main phonon 
lines below loq and the weak orbiton induced satellites at 
« 2uJo are in excellent agreement with experiment. The 
orbiton satellites are, just as in the localized limit (see 
Fig. Hp), due to the mixing of orbital and phonon modes. 
The value of the e-p coupling that is used in the phonon 
calculation, g/wQ — 0.35, corresponds to rather weak 
electron-phonon coupling, in accordance with Ref. p7[ , 
but in contrast to Ref. ||l^. If we were to use a larger 
value of g, the orbiton dispersion would become too 
small. A way to determine the e-p coupling regime exper- 
imentally is to check for additional satellites in the Ra- 
man spectrum at about w « 3luq, which for the coupling 
strength in the present calculation have low intensity, 
but would have large intensity if the system were in the 



the features around 150 meV in recent Raman-scattering 
experiments on LaMnOa as orbiton derived satellites in 
the phonon spectral function. These satellites should also 
be observable in other experiments that probe phonon 
dynamics, for instance in neutron scattering. 



strong coupling regime |12|. It is crucial that we took the 
phonon-dynamics into account, as a purely static Jahn- 
Teller distortion (as in Ref. |1^ , see also Fig. ||c) would 
lead to an orbital dispersion that is about a factor two 
too large. We interprete the Raman peaks around 150 
meV as orbiton satellites of the phonon peaks around 
80 meV, which would disappear without electron-phonon 
coupling. This is in contrast to Saitoh et al. [Q, where 
the peaks are due to different Raman scattering mecha- 
nisms and are assumed to be independent, which implies 
strong dependence of their intensity ratios on incident 
photon energy. 

Conclusions. We calculated the orbiton and phonon 
properties for a realistic model Hamiltonian for LaMnOa 
and compare the results with Raman-scattering data. 
We treat electron-electron correlations and lattice dy- 
namics on equal footing and have shown that the or- 
biton dispersion, which is mainly caused by correlation 
effects, is strongly reduced by the electron-phonon cou- 
pling. This coupling also mixes the orbiton and phonon 
modes and causes satellite structures in the orbiton and 
phonon spectral function. The elementary excitations of 
the system, in other words, are mixed modes with both 
orbital and phonon character. This leads us to interprete 



[1] M. Imada, A. Fujimori and Y. Tokura, Rev. Mod. Phys. 

70, 1039 (1998), Y. Tokura, N. Nagaosa, Science 288, 

462 (2000). 
[2] K. I. Kugel and D. I. Kliomskii, Sov. Phys. JETP 37, 

725 (1973); Sov. Pliys.-Usp. 25, 232 (1982). 
[3] H. A. Jalin and E. Teller, Proc. R. Soc. A161, 220 (1937). 
[4] H. Ichikawa et al, Physica B 281, 482 (2000), H. Sawada 

and K. Terakura, Phys. Rev. B. 58, 6831 (1998). 
[5] C. Castellani, C. R. Natoli and J. Ranninger, Phys. Rev. 

B. 18, 4945 (1978). 
[6] H.F. Pen et al, Phys. Rev. Lett. 78, 1323 (1997). 
[7] Y. Ren et al, Nature 396, 441 (1998). 
[8] J. B. Goodenough, Magnetism and Chemical Bond (In- 

terscience, New York, 1963). 
[9] H. Kuwahara et al, Science 270, 961 (1995). 
[10] S. Ishihara and J. Inoue, Phys. Rev. B 55, 8280 (1997). 
[11] J. van den Brink et al, Phys. Rev. B 59, 6795 (1999). 
[12] P.B. Allen and V. Perebeinos, Phys. Rev. Lett. 83, 4828 

(1999). 
[13] E. Saitoh et al, Nature 410, 180 (2001). 
[14] P.B. Allen and V. Perebeinos, Nature 410, 155 (2001). 
[15] L. F. Feiner, A. M. Oles, and J. Zaanen Phys. Rev. Lett. 

78, 2799-2802 (1997). 
[16] J. van den Brink et al, Phys. Rev. B 58, 10276 (1998). 
[17] J. van den Brink, P. Horsch and A.M. Oles Phys. Rev. 

Lett. 85, 5174 (2000), W.G. Yin, H.Q. Lin and CD. 

Gong, Phys. Rev. Lett. 87, 047204 (2001). 
[18] V. Perebeinos and P.B. Allen, Phys. Rev. Lett. 85, 5178 

(2000). 
[19] T. Hotta et al, Phys. Rev. B 60 R15009 (1999). 
[20] G. KhaliuUin and R. Kilian, Phys. Rev. B 61, 3494 

(2000). 
[21] Note that the 'mixed terms' in Horb, containing prod- 
ucts T^T^, are neglected. This is a reasonable approx- 
imation ]ll| , [17[ . This also holds for terms of the same 

symmetry in Hph- 
[22] J. Bala and A.M. Oles, Phys. Rev. B 62, R6085 (2000). 
[23] See, for instance, G.D. Mahan, Many-Particle Physics, 

Plenum, New York (1990). 
[24] J. van den Brink, to be published. 
[25] C. L. Kane, P. A. Lee and N. Read, Phys. Rev. B 39, 

6880 (1989), G. Martinez and P. Horsch, Phys. Rev. B 

44, 317 (1991). 
[26] A. Ramsak, P. Horsch and P. Fulde, Phys. Rev. B 46, 

14305 (1992). 
[27] P. Benedetti and R. Zeyher, Phys. Rev. B 59, 9923 

(1999). 



